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We propose an experimental scheme to create spin-orbit coupling in spin-3 Cr atoms using Raman 
processes. By employing the linear Zeeman effect and optical Stark shift, two spin states within 
the ground electronic manifold are selected, which results in a pseudo-spin- 1/2 model. We further 
study the ground state structures of a spin-orbit-coupled Cr condensate. We show that, in addition 
to the stripe structures induced by the spin-orbit coupling, the magnetic dipole-dipole interaction 
gives rise to the vortex phase, in which a spontaneous spin vortex is formed. 

PACS numbers: 37.10.Vz, 03.75.Mn 



Over the past few years, there has been rapidly grow- 
ing interest in engineering Abelian and non-Abelian arti- 
ficial gauge fields in ultracold atomic gases . Partic- 
ularly, the non-Abelian gauge field, or more specifically 
the spin-orbit (SO) coupling, is of fundamental impor- 
tance in many branches of physics. Fascinating examples 
include the quantum spin-Hall effect and the topologi- 
cal insulators in condensed matter physics Q. With the 
enormous tunability of the interaction and geometry, ul- 
tracold atomic gases may offer a tremendous opportunity 
for studying exotic quantum phenomena in many-body 
systems with SO couphng [sl-llq. 

In their pioneer experiments, the NIST group have re- 
alized the light- induced vector potentials [13, the syn- 
thetic magnetic fields ^Ts*], and the electric forces [l9| in 
ultracold Rb gases through Raman processes which 
differs from most dark-state based theoretical propos- 
als [2^ in that the linear Zeeman shift is compensated 
by the two-photon detuning. More remarkably, they also 
created a two-component SO-coupled condensate of Rb 
atoms and observed the phase transition from spatially 
mixed to separated states [2ll |. An important ingredient 
in this experiment is that the quadratic Zeeman shift is 
employed to separate two desired spin states from the re- 
maining one. Hence, this scheme is inapplicable to atoms 
without nuclear spin, such as certain isotopes of Cr and 
Dy, in which the quadratic Zeeman effect is absent. 

In this Letter, we propose an experimental scheme to 
create SO coupling in spin-3 ^^Cr atoms by selecting two 
internal states from the J = 3 ground electronic mani- 
fold. Similar to the NIST group's scheme, ours also re- 
lies on Raman processes. However, we utilize the optical 
Stark shift to compensate the linear Zeeman shift so that 
the lowest two levels are near degenerate and well sepa- 
rated from other levels, which leads to a pseudo spin- 1/2 
model. The proposed scheme has the advantages that 
only a moderate magnetic field strength is required and 
it also applies to atoms without nuclear spin. 

An interesting feature of the Cr atom is that it pos- 
sesses a large magnetic dipole moment, which makes the 
scalar Cr condensate an important platform for demon- 
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FIG. 1: (color online), (a) Scheme for creating SO coupling in 
a Cr atom. Two Raman beams, propagating along x + y and 
—X + y with frequency difference Aljl, are linearly polarized 
along z and x + y, respectively. A bias field Bo is applied 
along the negative z axis, which generates a Zeeman shift 
Ldz in the ground state manifold, (b) Level diagram for the 
Raman coupling within the \ J — 3) ground state manifold by 
utilizing the | J' ~ 2} excited state. 



strating the dipolar effects [22|. Moreover, when an 
atom's spin degree of freedom becomes available, mag- 
netic dipole-dipole interaction (MDDI) also couples the 
spin and orbital angular momenta, which is responsible 
for the Einstein-de Haas effects [l^, 01 , the spontaneous 
demagnetization |25|| of the Cr condensate, and the spon- 
taneous spin vortices [26l - l28j in spinor condensates. Un- 
fortunately, in spin-3 Cr condensates, contact interaction 
also contains spin-exchange terms which are much larger 
than the strength of the MDDI [2^ . Therefore, the spin 
vortex phases are yet to be observed. In the pseudo spin- 
1/2 Cr condensate, we show that only the MDDI con- 
tains spin-exchange terms and a spontaneous spin vortex 
is readily observable. 
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We consider a condensate of ^^Cr atoms subjected to 
a bias magnetic field Bg along the negative z-axis. The 
Zeeman shift within the ground state manifold is hujz — 
f/s/islBol with gs — 2 being the electron spin g-factor 
and hb the Bohr magneton. Here, the quadratic Zeeman 
shift is zero because of the absence of the nuclear spin. 
As shown in Fig. [1] atoms are illuminated by a pair of 
linearly polarized Raman beams which propagate along 
X + y and —x + y with frequencies ujl + Awl and ujl, 
respectively. The ground- (^5*3) to excited-state CP2) 



transitions are coupled by the Rabi frequencies flie'^^^ '^ 
and ri2e*'^2 '", where ki = kL{x + y) and k2 — kL{—x + y) 
are the wave vectors of the Raman beams with = 
v^tt/A and A being the wave length of the lasers. For 
simplicity, are assumed to be real. If the frequency of 
the lasers is far detuned from the ground- to excited-state 
transition, i.e., |rii.2/A| ^ 1 with A being the detuning, 
the excited states can be adiabatically eliminated to yield 
the atom-light interaction Hamiltonian 
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where Ac = ujz + Awl is the two-photon detuning, 
VLr = ~VLiVl2/A is the Rabi frequency for the Raman 
coupling, Ui^2 = — f^i2/A are the optical Stark shifts 
induced by the laser fields fii and f22, respectively, and 
T{t) = e^*^"^* and X{x) = e^**^^^ are introduced for 
short-hand notation. The physical significance of Eq. ([1]) 
can be readily understood '3Ci] by using the level diagram 
[Fig.IUb)]. 

From Hamiltonian ([T]), it is apparent that, under the 
conditions Ac -I- ?7i « and |J72|, \^r\ ^ |Ac|, the en- 
ergy levels mj — 3 and 2 can be separated from other 
levels due to the large Zeeman shift. These conditions 
can be satisfied by choosing ujz — /S. and assuming 
that IAwl/w^I <C 1 and |f22/^^i| ^ 1, which eventually 
leads to an effective two-level Hamiltonian: 
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for pseudo spin- up | t) = I"! J = 3) and -down | \) = 
1 2), where / is the identity matrix and a constant term, 
-(C/2 -I- Aujl/2)1, has been added to obtain Eq. We 
note that the atom-light interaction term in h can be 
intuitively treated as an effective magnetic field, 

Boff = h{gsHB)^^{2i^RCOs2kLX, —2QRsin2kLX, —Aujl)- 

Unlike the NIST group's scheme [Hj, here, an optical 
Stark shift — 57^/A is used to compensate the linear Zee- 
man shift, so that only the levels mj — 3 and 2 are 
Raman coupled near resonance (Awl ~ 0). 

To proceed further, let us focus on the motion of an 
atom along the x axis by freezing its y and z degrees of 
freedom. By applying a simple gauge transform [l2| , the 



single-particle Hamiltonian can be recast into 

/ ^2q2 \ ^ hAuji, 
K = [ + ElJ i + 2Kqaz + hQ.Ra^ 2"^*^^' 

where q = Px/f^ is the quasimomentum. El = 
^k\/ {2M) is the single-photon recoil energy, a-^.y.z are 
the Pauli matrices, and k = E^/kL is the SO coupling 
strength. Even though k is independent of Raman cou- 
pling strength, SO coupling strength is still tunable by 
varying the relative angle of the Raman beams It 
can be readily shown that, after dropping the constant 
El term, the eigenenergies of Eq. ([3]) are 
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in analogy to those in the spin-1 Rb condensate. In par- 
ticular, on the lower branch i?_(g), there exist two local 
minima at q± ~ zLkL\/i — /i^f2^/(4i?£) when hflR < 
2El and HAwl ^ El- The corresponding energies are 
E.{q±) ~ -El - h^Vty{AEL) ± hAujL/2. The states 
with quasimomenta fi.g_ and Hqj^ (labeled as | t') and 
I J,'), respectively) represent the dressed spin states in 
which atoms condense in the absence of the interactions. 

Here, we would like to discuss the experimental feasi- 
bility of our scheme. The transition wavelength from the 
ground to excited state is 429. Inm, which corresponds to 
a recoil energy EL/h — (27r)10kHz. Other laser param- 
eters can be set up as follows. Since the linear Zeeman 
shift in our proposal plays the role of the quadratic 
Zeeman shift in the NIST experiment [2l| . we may set 
hioz — 3.8El, which implies that the laser intensity 
l^lip — 3.8EL\A\/h is about the same order of magni- 
tude as that used in the experiment. To allow the Ra- 
man coupling D,R to vary from to El, which covers the 



3 



most interesting parameter region in the experiment, the 
maximum value of |f22| can be chosen as 0.26|f2i|. Con- 
sequently, the maximum value of U2 is less than 0.26i?L, 
which justifies the neglecting of U2 in Eq. ([T]). Finally, we 
point out that the SO coupling strength k in our scheme 
is 3.14 times larger than that in the Rb experiment due 
to the smaller mass and the shorter transition wavelength 
of Cr atom. 

Now we turn to study the many-body effect in a SO- 
coupled Cr condensate. To this end, we first write down 
the single-particle Hamiltonian, which, in the second 
quantized, takes the form 



where V{r) 



h + V{r)-^i *(r), 



(5) 
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7 z )/2 is an axially 



symmetric harmonic trap with uj± being the radial trap 
frequency and 7 the trap aspect ratio, ji is the chem- 
ical potential, and ^'(r) = [il;^{r),^i{r)]'^ is the field 
operator for the bare spin states. We note that Hq 
can also be expressed in terms of dressed spin states 
by using the transform ip^{r) ~ '0t'('") ~ ee^*'^'-^-!/'^' (r) 
and ipi{r) ~ —tpi'{r) + £e~^*'^-^^-!/'t'('')i where e ~ 
h^lfj./ {^El + HAlul) <C 1 in the weak Raman coupling 
limit h£tR/EL < 1. 

In terms of the bare spin states, the coUisional inter- 
action takes the form 
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where 54.6 = A-Kh^ajifi/M with 04 = 58 and ag = 
112 as being the s-wave scattering lengths for the coUi- 
sional channel with total spin angular momentum j = 4 
and 6, respectively [3l|. This result can be understood 
as follows. The collision between two mj — 3 atoms 
can happen only in the total spin j — 6 channel; con- 
sequently, it has a scattering length ag. For collisions 
between mj = 3 and 2 atoms, the projection of the to- 
tal spin along the z-axis is rrij = 5, which is conserved 
during collision. Therefore, this collision also happens in 
the J = 6 channel. But when two mj — 2 atoms collide 
with each other, both j — 6 and 4 channels will con- 
tribute. Apparently, the spin-up and -down states are 
immiscible. 

The MDDI for the pseudo spin-1/2 system can be de- 
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FIG. 2: Integrated densities (columns 1 and 3 for spin-up 
and -down, respectively) and phases of the condensate wave 
functions on the z = Q plane (columns 2 and 4 for spin-up 
and -down, respectively). From the first to the fourth rows, 
the frequency differences are H^ujl/El = 0.01, 0.04, 0.0875, 
and 0.1, respectively. 



here, ga = /io5s/i^/(47r) with /^o being the vacuum per- 
meability and jiB the Bohr magneton, e — (r — r')/|r — r'| 
is an unit vector, and we have adopted the notations 
V^Q = ■0Q(r) and ^'^ = ipai^') with a =t and 4-. The 
first three terms of H^^^ represent the intra- and inter- 
species dipolar interactions in a mixture of raj — 3 and 2 
atoms, and the last term is the exchange dipolar interac- 
tion. ' is of particular interest. It represents the SO 
coupling containing in the MDDI and does not conserve 
the atom number in the individual spin state. However, 
the total angular momentum is conserved by i?^ . 

In this spin-1/2 model, interactions have a much sim- 
pler form compared to those in the spin-3 system. In 
particular, here, only the MDDI contains spin-exchange 
terms. As will be shown, even though is much smaller 
than 34,6, the spin associated dipolar effect can be readily 
detected in pseudo spin-1/2 Cr condensates. 

We now investigate the ground state structures of the 
SO-coupled dipolar condensate using the mean-field the- 
ory. To this end, the field operators t\}a are replaced 
by the condensate wave function — {'tpa), which can 
be obtained by numerically minimizing the free energy 
functional V'i] — {Hq + Hs + Hd). Specifically, we 

consider a Cr condensate with iV = 10^ atoms. The pa- 
rameters for the trapping potential are chosen as a;_L — 
(27r)100Hz and 7 = 6, representing a three-dimensional 
pancake-shaped trap. Furthermore, the Rabi frequency 
for Raman coupling is fixed at Mln = — O.Oli^L. Since 



4 



a 

a. 



-4 




-8 -4 4 8 

FIG. 3: (color online). Vector plot of the transverse compo- 
nents of s(r) on the z — plane for HAlul ~ 0.0875El. The 
grayscale indicates the integrated density nf{x,y). 
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FIG. 4: (color online). Reduced atom number Na (a) and av- 
erage orbital angular momentum Lz (b) as functions of Aujl- 



£ ^ 1 is satisfied, we shall discuss only the ground state 
in terms of the bare spin states. 

In Fig. [2l we plot the integrated density, na{x,y) = 
J (iz|V'a(r)P, and the phases of the condensate wave 
functions for various Awl's. When Aujl is small (the 
first row in Fig. [5]) , the single-particle energies of the two 
pseudo spin states are nearly degenerate such that the 
ground state structure is mainly determined by the in- 
teractions. Apparently, both and Hj^ (in the pancake- 
shaped trap) favor the spin-down state. Therefore, | J, ) 
becomes dominantly populated, which we refer to as the 
polarized phase (PP). As shown in the fourth row of 
Fig. [21 the PP also occurs when the frequency difference 
Acji (or, equivalently, the z component of the effective 
magnetic field Boff) is sufficient large, under which | f) 
becomes dominantly occupied. 

The common feature of the PPs is that the wave func- 
tion of the highly populated state is structureless, as 
any structure developed in the high density spin state 
would cost too much kinetic energy. On the other hand, 
a striped structure forms in both the density and phase 
of the less populated spin state. The phase stripe can be 
intuitively understood as follows. To lower the energy, 
the pseudo spin density, s(r) — J^a/i i^a^aisfpp, has to be 
antiparallel to the local effectively magnetic field, which 
requires the relative phase of the condensate wave func- 
tion to take the form arg(-0^) — arg(-0^) ^ Tr+2kLX. Since 
the phase of the highly populated state is a constant, the 



phase of the other spin state is then periodically modu- 
lated along the x direction. The density stripe in the less 
populated spin state is caused by the immiscible nature 
of the two-component condensate. 

More remarkably, we observe a vortex phase (VP) for 
intermediate Awl values. As shown in the second and 
third rows of Fig. [2l a singly quantized vortex appears in 

" (2) 

the spin-up state due to the SO coupling in the term 
of the MDDI. In the VP, the atom numbers in spin-up 
and -down states become comparable [Fig. lUa)], which 
allows the spin s(r) of the atom to form significant trans- 
verse components. As shown in Fig. O since the MDDI 
is minimized with a head-to-tail spin configuration, the 
transverse components of s(r) are forced to form a spin 
vortex. Consequently, the wave function t/j^ develops a 2tt 
phase winding, representing a vortex state. The reason 
that the vortex state appears only on the spin-up com- 
ponent is due to the immiscibility of our two-component 
system, which results a density depletion at the center 
of V't- Therefore, forming a vortex in the spin- up state 
costs less kinetic energy. Moreover, in the VP, the phase 
stripes also appear in the low density regions of both spin 
states, which is the manifestation of the SO coupling in- 
duced by the light fields. 

To determine the phase boundaries, we plot the 
Alul dependences of the reduced atom number, = 
N^^ J drl'ipal'^, and the average orbital angular momen- 
tum, = I rfri/'^L^V'a, in Fig. H where = 
—ih{x-^ — y-§^) is the z-component of the orbital angu- 
lar momentum. As can be seen, the phase boundaries are 
marked by two critical Aujl values, Aw^ = 0.031i?L/S 
and Acj2* = OmSEL/h. For Aw^ < Awl < Aw£*, 
the condensate lies in the VP; otherwise, it is in the PP. 
Within the VP, atom numbers and orbital angular mo- 
mentum change dramatically as one varies Auj^. 

In conclusion, we have proposed an experimental 
scheme to generate SO coupling in spin-3 Cr condensates 
via Raman processes. Optical Stark shift is employed to 
selecting two spin states from the atom's ground elec- 
tronic manifold. The proposed scheme should be read- 
ily realizable experimentally. Subsequently, the ground- 
state structures of a SO-coupled Cr condensate have been 
investigated. We show that the interplay between the 
light fields and the MDDI gives rise to the polarized and 
vortex phases. In particular, a spontaneous spin vortex 
is formed in VP. The spin vortex state can be experimen- 
tally identified if it is observed that the spin-up conden- 
sate is a vortex state and the spin-down condensate is a 
vortex-free one. Finally, we point out that our scheme 
should also apply to the Dy atom ^] , which has an even 
larger dipole moment. 
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